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Abstract 

The isotopic shifts are determined to high accuracy for a number of hght two-electron Li+, Be^+, 
B'^^ and C^^ ions. The field components of these isotopic shifts have been determined with the 
use of the exact Racah-Rosental-Breit formula. 
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The dependence of the total energies of bound states of different isotopes of the same 
chemical element upon the 'isotope numbers' is called the isotopic shift. By the 'isotope 
numbers' we mean the mass of the nucleus and its electric charge density distribution (or 
proton density distribution). In general, the isotopic shift can be observed in arbitrary, in 
principle, atomic system which has a central, heavy nucleus with the finite mass, positive 
electric charge and finite radius. In general, the isotopic shift AE of the bound state 
level with the total energy E has the two following components: the first component which 
depends upon the mass of the nucleus and the second component which mainly depends upon 
the electric charge density distribution in the atomic nucleus. In turn, the first component is 
represented as the sum of the normal and specific components. Each of these two components 
is proportional to the factor where me is the mass of the electron, while M is the nuclear 
mass. For few-electron (A^— electron) atoms and ions the exact formula for the isotopic shift 
AE takes the form 

N 2 N 

ae = ae^ + ae, = ^{J:^) + ^{ E Pp^) (1) 

where the notation {X) designates the expectation value of the operator X. For the two- 
electron (or helium-like) atoms and ions in atomic units h = l,me = l,e = 1 the last formula 
is reduced to the form 

AE = AE„. + ^E, = ^{pI) + ^ (pi ■ P2) (2) 

As follows from this formula the expression for the normal and specific components of the 
isotopic shift are 

A^n = ^(p?) and Ai?, = ^(pi-p2) , (3) 

respectively. As follows from Eq.Q to determine the normal and specific components of 
the isotopic shift in two-electron atom/ion one needs to obtain the expectation values of the 
pI and pi ■ p2 operators. Everywhere in this study we assume that the wave functions of 
the two-electron atom/ion are properly symmetrized upon spin-spatial permutations of the 
two electrons and, therefore, the corresponding single-electron expectation values are always 
equal to each other, e.g., (p^) = (P2)- 

In actual two-electron atomic systems, i.e. in the systems with the finite nuclear mass 
M, one can use the condition which follows from the conservation of the total momentum 



Pat = Pi + P2, where Pat is the momentum of the nucleus, while pi and p2 are the electron 
momenta. From here one finds: 

^(P^) = (P?) + (Pi-P2) (4) 
and, therefore, from Eqs.([2]) and we have 

(6) 

i.e. the mass dependent component of the isotopic shift is the expectation value of the 
kinetic energy of the atomic nucleus with the large (but finite!) mass. In many books and 
textbooks the formula, Eq.([5]), is considered as the original (or fundamental) expression, 
while Eq.([T]) is derived from this formula. 

As mentioned above the field component of the isotopic shift explicitly depends upon the 
nuclear size (radius) R and proton distribution density in the nucleus. It is clear that this 
component also depends upon the mass of atomic nucleus M, since the nuclear matter is 
a saturated matter (in contrast with the Coulomb matter). Therefore, the nuclear radius 
R is uniformly related to the number of nucleons A in the nucleus: R = Vq ■ A^, where 
the 'constant' radius tq ~ 1.17 — 1.25 ■ 10^^'^ cm = 1.17 - 1.25 fm (fermi), where 1 fm = 
1 ■ 10^^^ cm. Briefiy, this means that the field component of the total isotopic shift is also 
a function of the nuclear mass M, since A ^ where is the proton mass. The more 



accurate formula for the nuclear mass as the function of A, Z{= Np), where Z is the electric 

charge of the nucleus = number of protons Np) and (number of neutrons) is given by 

the Weizacker formula. This formula is discussed in the Appendix. 

Let us present here the well known formula obtained by Racah, Rosental and Breit for 

the field shift (see, e.g., jlj and references therein) 

Anal ('^QRx^b 6R 

^^ = ^- [r(2fe + i)p -'^(^^-i^) ■^^■^^(r-)) (6) 

where Q is the nuclear charge, R is the nuclear radius and b = y/1 — a'^Q'^, where a = 
^ !^ ^ is the dimensionless constant which is the small parameter in QED. In Eq.(l3]) 
the notation T{x) stands for the Euler's gamma-function, while the factor B{b) is directly 
related to the proton density distribution in the atomic nucleus. By assuming the uniform 
distribution of the proton density over the volume of the nucleus one finds the following 
expression for the factor B{b) 



For light nuclei with Q < 6 we have b ^ 1 and i? ~ | . The formula, Eq. (|6]) , was used in many 
theoretical papers for numerical evaluations of the filed component of the isotopic shift, or 
field shift for short. In some works, however, this formula was written with a number of 
'obvious simplifications'. Many of such 'simplifications' are based on the fact that for light 
nuclei the numerical value of the factor b is close to unity. Furthermore, in some papers the 
factor b was mistakenly called and considered as the Lorentz factor, while the actual Lorentz 
factor 7 is its inverse value 7 = = "^y=^^- Such a factor 7 always exceeds unity. 

In this study we evaluate the isotopic shift for a number of the ground 1^S{L = 0)— states 
in the light two-electron ions by using the exact (not approximate!) formula, Eq.([n]). Such 
evaluations allow one to evaluate the numerical errors which arise from the use of approx- 
imate expressions. Our main interest in this study is related with the field component of 
the isotopic shift. As follows from Eq.([6]) to evaluate the field component of the field shift 
one needs to determine to very high accuracy the expectation value of the electron-nuclear 
delta-function (or the {6{reN)) value). In this study for each light element/atom all isotopic 
shifts are determined in respect to the model (or idealized) isotope which has the infinitely 
heavy nucleus with zero spatial radius (zero charged radius). In this case the ratio ^ in 
Eq.(l6]) equals unity, and this equation can be re- written to the form (in atomic units) 

^M=^ -Q •« ■[r(2fe+l)P(26+l)(26 + 3)-y -^^('^^^^ 

where re = a^ao ~ 2.817940 fm (1 fm = Ic- 10~^^ cm)) is the classical radius of the electron. 
For light atomic nuclei the dimensionless factor ^ in the last formula is close to unity. Also, 
in our calculations we have used the following numerical values for the physical constants: 
a = 7.297352569 ■ lO'^ and ao = 5.291772109 ■ 10"^ cm. The formula, Eq.(^, has been used 
in all calculations performed for this study. As follows from Eq.([H]) to determine the field 
component of the isotopic shift one needs to know the expectation value of the electron- 
nuclear delta-function 6{reN) and numerical value of the nuclear radius R. The expectation 
value of S^Tsn) can be found from the results of highly accurate atomic computations, while 
the nuclear radius must be taken from nuclear experiments (see, e.g., Jsl). 

In this study we consider a few light two-electron (or He-like) ions: Li+, Be^^, and 
C^"*". The total energies and some other properties (or expectation values) of the corre- 
sponding model ions (with the infinite nuclear mass and zero spatial radii) can be found in 
Tables I - IV. The convergence of the results upon the total number of basis function 
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used in calculations is also shown in Tables I - IV. For these Tables we have computed the 
expectation values of the following operators: 5(reAr), i^eN, Pi, Pi ■ P2 and p^. The notation 
i^eN designates the electron-nuclear cusp value which is defined by the following equation 

Formally, it is an avaraged velocity of the electron at the nucleus. In the general case, this 

expectation value gives the relative velocity of the two particles {i and j) at the (^j)— collision 

point. For pure Coulomb systems, e.g., for the atom with the nuclear charge Q, such a 

velocity is known from the corresponding classical problem. In particular, for the atom with 

the nuclear charge Q, the electron-nuclear cusp u^n must be equal (in atomic units) 

2 rUeMN . . 

T^eN = -Qe -—- = -Q 10) 

rUe + Mn 

since = 00. 

The coincidence of the computed z/gjv value with the nuclear charge Q (or -Q) indicates the 
quality of the variational (atomic) wave functions around the nucleus and overall accuracy 
of the computed expectation value of the electron-nuclear delta-function. In all calculations 
performed for this study we have used our exponential variational expansion in relative 
coordinates described in detail in our earlier papers, see, e.g., 

a, a 

. Here we do not want to 

repeat these descriptions of our variational expansion. Note also that in our calculations we 
have determined many dozens of different expectation values, including singular expectation 
values. However, in Tables I - IV only to a few expectation values are presented. All these 
values are needed for numerical evaluation of the isotopic shifts. 

Table V contains the numerical values of the field components of isotopic shifts (in a.u.) 
derermined with the use of the formula Eq.Q. In this Table we also present the numerical 
values of the following factors from that formula: R (the actual nuclear radius), b,X = 
Af'^^nQ'^''^^ ■ a^'' ■ [r(2b+i)^^l2b+i)(2b+3) ^ ~ (r^) • "^^^ expectation values of the electron- 
nuclear delta-functions were taken from Tables I - IV. To evaluate the Euler's gamma- 
function T(x) we have used the approximate 7-term formula derived by Lanczos [7|. The 
overall accuracy of this formula is ~ 1 ■ 10"^*' — 2 ■ 10^^'^. 

Appendix I. 

The formula which provides the uniform relation between the nuclear mass M and total 
number of nucleons A, nuclear charge Z (= number of protons Np) and number of neutrons 



Nn in the nucleus was derived in 



i J. 

a 



937 by Bethe, Weizacker and others. Now, this formula is 



know as the Weizacker formula |8|], or Bethe- Weizacker formula. This five-term formula for 
the nuclear binding energy Ei, was produced 75 years ago and since then its general structure 
has never been changed. First, note that the mass formula for an arbitrary nucleus is written 
in the form 



M 



^ + iV(-)-A] (11) 



where M is the nuclear mass of the nucleus with A nucleons, Z protons and neutrons 
A = Z + N). Also in this formula Eb is the binding energy of the nucleus, c is the speed of 
light in vacuum, while nip and m„ are the masses of the proton and neutron, respectively. 
The factors rUpC^ = 938.272910 MeV and m„c2 = 939.565378 MeV. The advantage of 
the formula, Eq. ffTT]) . is obvious, since it contains only dimensionless ratios and two integer 
numbers {Z and A^). For instance, if we chose in Eq.f llip rup = 1836.152701me, then M will 
be given in rrie (or in atomic units if rrie = 1). This is very convenient for highly accurate 
computations of different few-electron ions. 

The parameter Ei, in Eq.( lTT|) is called the binding energy of the nucleus. The explicit 
expression for the nuclear binding energy Eh is written as the following sum (the Weizacker 
formula) : 

2 Z"^ ( N — Z)'^ 

Eb = avA-asA^ -ac—-aA- —^ + S{A,Z) (12) 

As A 

where the five terms in the right-hand side of this equation are called the volume term, 
surface term. Coulomb term, assymetry term and pairing term, respectively. The pairing 
term 6{A, Z) equals zero, if A is odd. If A is even and both Z and are even, then 
6{A, Z) = The Weizacker formula is relatively accurate for regular nuclei (i.e. for nuclei 
which are not far from the stability region). In reality, such an accuracy directly depends 
upon the numerical values of parameters ay, as, ac, a a and Op in Eq. flT^ . In our calculations 
we have used the following values of these parameters: ay = 15.8 MeV , as = 18.3 MeV, 
ac = 0.714 MeV, aA = 23.2 MeV and ttp = 12.0 MeV. The Weizacker formula with these 
coefficients is sufficiently accurate for all light nuclei which are located in the stability region. 

Appendix II. 

In general, the bound state in any few-electron atom arises in the result of proper balance 
between kinetic energy Tg of atomic electrons, electron-nuclear (Coulomb) attraction Uen and 



electron-electron (Coulomb) repulsion Uee- The total energy E of this bound state is written 
as the following sum: 

E = T, + Une + Uee = ^ (^ne + U^e) (13) 

where in the right-hand side we have used the virial theorem for the Coulomb potential 
(see, e..g, 9|) Tg = —^{Une + Uee)- To solve the problem we need the explicit expression 
for the correlation energy Uee in terms of the electron- nuclear attraction energy Um- Such 
an analytical relation can be found in a few restricted cases, e.g., in the Thomas- Fermi 
method one finds: Uee = —\Une and E = ^Une- For highly accurate method the situation 
is significantly more complicated, but qualitatively the actual relations between Tg, Une and 
Uee are similar to the relations obeyed in the Thomas-Fermi model. 

First, note that the Schrodinger equation for an arbitrary non-relativistic A^— electron 
atom/ion takes the form 

%2 N TV -1 1 

^Ep?-Qe^E7 + e^ E ^]^ = E^ (14) 

^'"■e 1=1 1=1 ' t i>j=l 

where Q is the electric charge of the nucleus and the mass of the central nucleus is assumed 
to be finite. Here we shall not assume that N = Q. In Eq. ffT^ the terms in the first 
sum (~ p^) represent the kinetic energies of electrons, while other terms correspond to the 
electron-nuclear Coulomb attraction and/or electron-electron Coulomb repulsion. From this 
equation one finds the following relation for the expectation values (in atomic units): 

[f(p,^)-W^ + M(±)].^ (15) 

^ 'en ^ ' ee 

where the subscript e means the electron, while the subscript n stands for the nucleus. 

Let us assume that we know the highly accurate solution of the Schrodinger equation, 
i.e. the wave function in Eq.([I3j) is known. In this case the following condition must be 
obeyed for these three expectation values: 

'en ^ ' ee 

which is known as the virial theorem (see, e.g., [9i]). Now, after a few steps of simple 
transformations one finds the following formula for the total energy E (in atomic units) 

E = -In.Q.{—) + \n{N-1).{-) 

^ '"en ^ '"ee 

= -^N ■ Q ■ MN, Q) + -^N{N - 1) ■ /2(iV, Q) (17) 



where fi{N,Q) and f2{N,Q) are the functions of the number of electrons and electric 
charge of the nucleus Q. For neutral atoms we have Q = N, while for positively charged ions 
Q > N. For instance for the two-electron C"^^ ion (i.e. for = 2 and Q = 6) we found that 
(^) ^ 5.687615325399107305988929 a.n. and (^) ^ 3.438890700992227050791016158 a.?/, 
(results from calculations with 3700 basis functions). By using these expectation values one 
finds the total energy E from Table IV. In the general case, (^) = A(A^, Q){-^), where the 
paramter A < 1. The factor X{N,Q) can be approximated to very high accuracy with the 
use of highly accurate results obtained for different N and Q. 

The formula, Eq. f[T7|) . is also applied in the case of = 1 (hydrogen-like atoms and ions). 
In this case the electron-electron repulsion is not defined, but it is multiplied by a factor 
which equals zero. As follows from Eq. f[T7|) . if we can determine the two unknown functions 
fi{N,Q) and f2{N,Q) in Eq. ffTTl) . then we can predict the total energy of an arbitrary 
A^— electron atom/ion to very high accuracy. It can be done, e.g., by using the results of 
highly accurate computations of different few-electron atoms and ions (with different Q and 
A^). If we know the results of highly accurate computations for a large number of light 
atoms/ions, then we can reconstruct the two functions fi{N, Q) and f2{N, Q) in Eq.( |T71) to 
very high accuracy. This procedure works well for mixtures of light elements at different 
temperatures and allows one to predict the properties of high-temperature plasmas which 
consists of such elements. 
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TABLE I: The total energies E and expectation values of the electron- nuclear delta-function SgN, 
electron-nuclear cusp z^eat and some other operators for the two-electron lithium ion Li+ (in atomic 
units). K is the total number of basis functions used. 



K 



E{U+) 



{Scn) 



I'eN 



3500 -7.279913 412669 305964 91708 
3700 -7.279913 412669 305964 91743 
3840 -7.279913 412669 305964 91766 
4000 -7.279913 412669 305964 91785 



6.8520 094343 431 
6.8520 094343 456 
6.8520 094343 460 
6.8520 094343 462 



-3.0000 00000 158 
-3.0000 00000 125 
-2.9999 99999 918 
-2.9999 99999 901 



K 



(Pl • P2) 



HPn) 



3500 3.63995 670633 465298 240 

3700 3.63995 670633 465298 241 

3840 3.63995 670633 465298 241 

4000 3.63995 670633 465298 242 



0.288975 786393 989535 661 
0.288975 786393 989535 661 

0.288975 786393 989535 662 
0.288975 786393 989535 662 



7.56888 919906 329532 141 
7.56888 919906 329532 143 

7.56888 919906 329532 144 
7.56888 919906 329532 145 



TABLE II: The total energies E and expectation values of the electron-nuclear delta-function SeN, 
electron-nuclear cusp fgAr and some other operators for the two-electron berillium ion Be^"^ (in 
atomic units). K is the total number of basis functions used. 



K 



E{Be'^+) 



I'eN 



3500 -13.65556 623842 358670 20757 
3700 -13.65556 623842 358670 20767 
3840 -13.65556 623842 358670 20772 
4000 -13.65556 623842 358670 20777 



17.1981 72544 645 
17.1981 72544 640 
17.1981 72544 638 
17.1981 72544 635 



-3.9999 99999 962 
-3.9999 99999 921 
-4.0000 00000 125 
-4.0000 00000 148 



K 



(Pl • P2) 



3500 6.82778 311921 179335 084 

3700 6.82778 311921 179335 086 

3840 6.82778 311921 179335 089 

4000 6.82778 311921 179335 091 



0.420520 303439 441862 Oil 
0.420520 303439 441862 010 
0.420520 303439 441862 009 
0.420520 303439 441862 009 



14.07608 654186 302856 368 
14.07608 654186 302856 369 
14.07608 654186 302856 370 
14.07608 654186 302856 371 
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TABLE III: The total energies E and expectation values of the electron-nuclear delta-function SeN, 
electron-nuclear cusp UeN and some other operators for the two-electron boron ion B^"*" (in atomic 
units). K is the total number of basis functions used. 



K 



3500 -22.03097 1580242 781541 65339 
3700 -22.03097 1580242 781541 65376 

3840 -22.03097 1580242 781541 65394 
4000 -22.03097 1580242 781541 65418 



34.758 743660 955 
34.758 743660 965 

34.758 743660 947 
34.758 743660 935 



-5.0000 0000 319 
-5.0000 0000 235 
-5.0000 0000 107 
-5.0000 0000 119 



K 



{Pi ■ P-i) 



3500 11.01548 579012 139077 100 

3700 11.01548 579012 139077 089 

3840 11.01548 579012 139077 086 

4000 11.01548 579012 139077 083 



0.552752 631642 101467 789 
0.552752 631642 101467 734 
0.552752 631642 101467 715 
0.552752 631642 101467 701 



22.58372 421188 488300 979 
22.58372 421188 488300 952 
22.58372 421188 488300 942 
22.58372 421188 488300 938 
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TABLE IV: The total energies E and expectation values of the electron- nuclear delta-function SgN, 
electron-nuclear cusp I'eN and some other operators for the two-electron carbon ion C^"*" (in atomic 
units). K is the total number of basis functions used. 



K 



^(C4+) 



3500 -32.40624 660189 853031 05527 
3700 -32.40624 660189 853031 05535 

3840 -32.40624 660189 853031 05539 
4000 -32.40624 660189 853031 05542 



61.443 578056 445 
61.443 578056 514 

61.443 578056 537 
61.443 578056 543 



-5.9999 99998 765 
-5.9999 99999 871 

-6.0000 00000 048 
-6.0000 00000 037 



K 



(pi • P2) 



3500 16.20312 330094 926515 523 

3700 16.20312 330094 926515 524 

3840 16.20312 330094 926515 525 

4000 16.20312 330094 926515 525 



0.685334 822135 598924 527 
0.685334 822135 598924 535 
0.685334 822135 598924 536 
0.685334 822135 598924 537 



33.09158 142403 412923 500 
33.09158 142403 412923 502 
33.09158 142403 412923 502 
33.09158 142403 412923 503 



TABLE V: The nuclear radius R {fm), parameter b, factors X and Y (see the main text) and field 
components of the total isotopic shift AE^^ (all values are in atomic units) for each isotope. 



is()l()i")(' 


Q 


R b 


X 


Y ARf" 




6Li 


3 


2.5385 0.99976034018621 


4.297056149289-10-^ 


0.81154478 2.389469390-10" 


-7 


^Li 


3 


2.4312 0.99976034018621 


4.297056149289-10"^ 


0.74440369 2.191782729-10" 


-7 


9Be 


4 


2.5180 0.99957389838248 


5.749782211793-10"^ 


0.79852685 7.896292299-10" 


-7 




5 


2.4278 0.99933413638122 


7.219245621776-10-^ 


0.74241787 1.862963331-10" 


-6 


iiB 


5 


2.4059 0.99933413638122 


7.219245621776-10-^ 


0.72909310 1.829527229-10" 


-6 


12c 


6 


2.4073 0.99904101579314 


8.7092766851788-10^8 


0.76867983 4.113429628-10" 


-6 


13c 


6 


2.4614 0.99904101579314 


8.7092766851788-10-8 


0.76315629 4.083871553-10" 


-6 


14c 


6 


2.5037 0.99904101579314 


8.7092766851788-10-8 


0.78958608 4.225305034-10" 


-6 



12 



